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ABSTRACT 

Fluctuations in the redshifted 21 centimeter emission from neutral hydrogen probe 
the epoch of reionization. We examine the observability of this signal and the im- 
pact of extragalactic foreground radio sources (both extended and point-like). We use 
cosmological simulations to predict the angular correlation functions of intensity fluc- 
tuations due to unresolved radio galaxies, cluster radio halos and relics and free-free 
emission from the interstellar and intcrgalactic medium at the frequencies and angular 
scales relevant for the proposed 21cm tomography. In accord with previous findings, 
the brightness temperature fluctuations due to foreground sources are much larger 
than those from the primary 21cm signal at all scales. In particular, diffuse cluster 
radio emission, which has been previously neglected, provides the most significant 
foreground contamination. However, we show that the contribution to the angular 
fluctuations at scales 0^1 arcmin is dominated by the spatial clustering of bright 
foreground sources. This excess can be removed if sources above flux levels 5* > 0.1 
mjy (out to redshifts of z ~ 1 and z 2 for diffuse and point sources respectively) 
are detected and removed. Hence, efficient source removal may be sufficient to allow 
the detection of angular fluctuations in the 21cm emission free of extragalactic fore- 
grounds aX 9 > 1 arcmin. In addition, the removal of sources above 5 = 0.1 mJy also 
reduces the foreground fluctuations to roughly the same level as the 21cm signal at 
scales 6 1 arcmin. This should allow the substraction of the foreground components 
in frequency space, making it possible to observe in detail the topology and history of 
reionization. 

Key words: cosmology:theory - intcrgalactic medium - diffuse radiation - galaxies: 
active - radio continuum:general 



1 INTRODUCTION 

In the past few years a quantitative study of the high- 
redshift intergalactic medium (IGM) and its reionization 
history has finally been made possible by the discovery of 
quasars at a > 5.8 (e.g. Fan et al. 2001, 2003). In particular, 
the detection of a Gunn-Peterson trough (Gunn & Peterson 
1965) in the Keck (Becker et al. 2001) and VLT (Penter- 
icci et al. 2002) spectra of the Sloan Digital Sky Survey 
quasar SDSS 1030-0524 at 2 = 6.28 and in the Keck spec- 
trum of SDSS 1148-^5251 at 2 = 6.37 (White et al. 2003), 
has been interpreted as the signature of the trailing edge 
of the cosmic reionization epoch. The recent analysis of the 
first year of data from the Wilkinson Microwave Anisotropy 
Probe (WMAP) satellite on the temperature and polar- 
ization anisotropics of the cosmic microwave background 
(CMB), infers a mean optical depth to Thomson scatter- 
ing Te ~ 0.17, suggesting that the universe was reionized at 



higher redshift (Kogut et al. 2003; Spergel et al. 2003). Phys- 
ically, the CMB and the Gunn-Peterson trough probe two 
different stages of reionization, the former being sensitive to 
the initial phase, when free electrons appear, the latter to 
the residual neutral hydrogen in the latest stages of reion- 
ization. None of the two methods though is able to constrain 
the exact ionization level or the details of the reionization 
history. For this reason, an alternative way to probe the 
high-redshift IGM is required. 

An optimal experiment for probing the various stages 
of reionization is the proposed 21cm tomography. It has in- 
deed been shown that neutral hydrogen in the intergalactic 
medium (IGM) and gravitationally collapsed systems should 
be directly detectable in emission or absorption against the 
cosmic microwave background radiation (CMB) at frequen- 
cies corresponding to the 21cm line (e.g. Field 1958, 1959; 
Scott & Rees 1990; Kumar et al.l995). In principle it will 
possible to carry out such an experiment with planned high 
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sensitivity radio telescopes such as the PrimevAl Structure 
Telescope (PAST) the Square Kilometer Array (SKA) t 
and the LOw Frequency ARray (LOFAR) ^ . The 21cm spec- 
tral features will display redshift dependent angular struc- 
ture due to evolving inhomogeneities in the gas density field, 
hydrogen ionized fraction, and spin temperature. Several dif- 
ferent signatures have been investigated in the recent litera- 
ture: the fluctuations in the 21cm line emission induced both 
by the inhomogeneities in the gas density and in the ionized 
hydrogen fraction (Madau, Meiksin & Rees 1997; Tozzi et 
al. 2000; Ciardi & Madau 2003, hereafter CM; Furlanctto, 
Sokasian & tieriiquist 2004) and by 'minihalos' with virial 
temperatures below 10** K (Iliev et al. 2002, 2003); the global 
feature ('reionization step') in the continuum spectrum of 
the radio sky that may mark the abrupt overlapping phase 
of individual intergalactic HII regions (Shaver et al. 1999); 
and the 21cm narrow lines generated in absorption against 
very high redshift radio sources by the neutral IGM (Car- 
illi, Gnedin & Owen 2002) and by intervening minihalos and 
protogalactic disks (Purlanetto & Loeb 2002). 

While the 21cm tomography proposes to map the topol- 
ogy of the reionization process and constrain the nature of 
the ionizing sources, it remains a challenging project due 
to foreground contamination from unresolved extragalac- 
tic radio sources (Di Matteo et al. 2002), free-free emis- 
sion from the same halos that roionize the universe (Oh 
& Mack 2003) and the Galactic free-free and synchrotron 
emission (Shaver et al. 1999). Because the proposed exper- 
iments will be carried out in frequency space as well as an- 
gle, recent work has discussed the possibility of removing 
the foreground power spectrum components by comparing 
maps closely spaced in frequency (Zaldarriaga, Furlanetto & 
Hernquist 2004; Gnedin & Shaver 2004). The proposed sub- 
straction is nonetheless demanding as the foreground fluc- 
tuations overwhelm the primary 21cm signal by up to throe 
orders of magnitude. In particular, it will require knowing in 
detail the spectral behavior and possible spectral variations 
of each of the contaminants and their associated angular 
power spectra. 

In this paper we employ computer simulations of a 
ACDM universe to evaluate the foreground brightness tem- 
perature fluctuations due to extragalactic sources that con- 
taminate the 21cm tomography. In particular, we model the 
free- free emission from ionizing sources self-consistently with 
a viable reionization scenario and with the associated 21cm 
IGM emission (as described in CM). Using the same simula- 
tions we adopt simple but physically motivated prescriptions 
to model the radio galaxy population; the model success- 
fully reproduces both the observed luminosity function (at 
1.4GHz) and the two point correlation function of the ob- 
served radio sources. Finally, we use a separate simulation 
(of the same ACDM universe) that, in addition to dark mat- 
ter and baryonic gas, follows the evolution of magnetic fields 
and cosmic rays, to estimate the foreground signals due to 
extended radio sources such as cluster radio relics and radio 
halos, which had not been considered in previous work. 
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We shall show that the detailed information on the spa- 
tial and redshift distribution of the contaminant sources is 
crucial for determining the correct brightness temperature 
fluctuations due to the unresolved extragalactic foregrounds 
on the angular scales where the primary 21cm signal is ex- 
pected to peak. By modeling the foregrounds within the 
cosmological simulations we are able to study and map the 
variation of the angular clustering signal as a function of 
the flux cut above which foreground sources will be detected 
(and hence will not contribute to the unresolved fluctuating 
signal). In particular, we will show that the angular power 
spectra of unresolved foreground contaminants become sig- 
nificantly suppressed at scales > 1 arcmin when sources 
above fiux levels of a fraction of mjy are removed. 

The structure of the paper is as follows. In Section 2 
we describe the numerical simulations adopted to model the 
physical processes producing the extragalactic backgrounds. 
In Section 3 we briefiy outline the origin of the 21cm emis- 
sion line from the diffuse IGM used in CM, while in Section 4 
we show the total comoving emissivities and spatial corre- 
lation functions for the extragalactic foregrounds modeled 
within the simulations. In Sections 5, and 6 we calculate the 
contribution to the brightness temperature fluctuations to 
to the clustering of the different foregrounds and finally, in 
Section 7 we give our conclusions. 



2 NUMERICAL SIMULATIONS 

In this Section we briefly describe the numerical simulations 
adopted to model the various emission processes considered 
in this paper and refer to Ciardi, Stoehr & White (2003, 
hereafter CSW), Ciardi, Ferrara & White, (2003, hereafter 
CWF) and Miniati (2002) for further details. 

All simulations assume a ACDM "concordance" cos- 
mology with nm=0.3, Oa=0.7, h =0.7, ^6=0. 04, n=l and 
(78=0.9. These parameters are within the WMAP experi- 
mental error bars (Spergel et al. 2003). 

The simulations used to model the foregrounds pro- 
duced by radio galaxies and free-free emission from HII 
regions are based on the N-body code GADGET (Springel, 
Yoshida & White 2001). A "re-simulation" technique 
(Springel et al. 2001, hereafter SWTK) was employed in 
order to follow at high resolution the dark matter distribu- 
tion within an approximately spherical region with a diam- 
eter of about 50h~^ Mpc within a cosmological volume of 
479 Mpc on a side (Yoshida, Sheth & Diaferio 2001). 
The location and mass of dark matter halos was deter- 
mined with a friends-of-friends algorithm. Gravitationally 
bound substructures were identified within the halos with 
the algorithm SUBFIND (SWTK) and were used to build the 
merging tree for halos and subhalos following the prescrip- 
tion of SWTK. The smallest resolved halos have masses of 
M ~ 10^ Mq (the particle mass is Mp = 1.7 x W^h'^ Mq). 
The galaxy population was modeled via the semi-analytic 
technique of Kauffmann et al. (1999) and implemented in 
the way described by SWTK. We obtain in this way, a cat- 
alogue of galaxies for each of the simulation outputs, con- 
taining for each galaxy, among other quantities, its position, 
mciss and star formation rate. We used the galaxies within 
the 50h~^ Mpc region, coupled with the prescriptions de- 
scribed in Section 4, to model the redshift distribution and 



© 2004 RAS, MNRAS 000, 1-13 



21cm emission and foregrounds 3 



luminosity function of radio galaxy population and the free- 
free emission from the interstellar medium (ISM). 

Within the high resolution spherical subregion, a cube 
of comoving side L — 20h^^ Mpc was used to study the 
details of the reionization process, using the radiative trans- 
fer code CRASH (Ciardi et al. 2001; Maselh, Ferrara & Cia- 
rdi 2003) to model the propagation of ionizing photons into 
the IGM. Several sets of radiative transfer simulations were 
run in CSW and CFW, with different choices for the galaxy 
emission properties. The ones of used here axe those labeled 
S5 ('late' reionization case) and L20 ('early' reionization 
case) which employ an omission spectrum typical of Pop 111 
stars, a Salpeter Initial Mass Function (IMF) and an escape 
fraction of ionizing photons /esc = 5% and a Larson IMF 
with /esc = 20% respectively. For details and discussion on 
the choice of parameters we refer to CSW and CFWS. These 
simulations provide the redshift evolution and the spatial 
distribution of ionized IGM and have been used to model 
the free-free emission from the IGM. 

For the diffuse radio emission associated with cluster 
radio halos and relics we have used a different simulation of 
the same cosmological model, in which in addition to dark 
matter and baryonic gas, magnetic fields and cosmic-rays 
(CR) were also followed. In the simulation, the cosmological 
code of Ryu et al. (1993) is used in combination with the 
cosmic-ray code COSMOCR fully described in Miniati (2001, 
2002). Basically, during the simulation large scale structure 
shocks are identified and their jump conditions computed 
(Miniati et al. 2000). A fraction of the particles crossing the 
shock is assumed to be injected in the diffusive shock acceler- 
ation mechanism, and thereafter accelerated to a power-law 
spectrum in momentum space. The log-slope of the distribu- 
tion is determined by the shock Mach number (M) in accord 
with the test particle limit (e.g. Bell 1978). The fraction of 
injected CR protons is computed according to a variant of 
the thermal leakage prescription (e.g. Kang & Jones 1995); 
for practical purposes it amounts to about several xlO~®. 
The fraction of injected CR electrons is assumed to be 100 
times smaller than that (to be roughly consistent with mea- 
surements of this parameter from supernova remnants ob- 
servations and with EGRET 7-ray upper limits at 100 MeV 
for nearby clusters of galaxies; Miniati 2002). Both popula- 
tions of shock accelerated CRs (primaries) are then passively 
advected with the baryonic flow. In addition, we follow the 
evolution of their momentum spectrum by accounting for all 
relevant mechanisms of energy loss, such as Coulomb col- 
lisions, adiabatic, bremsstrahlung, synchrotron and inverse 
Compton losses for the electrons; and Coulomb and inelastic 
p-p collisions for the CR protons. As a byproduct of inelastic 
p-p collisions secondary are also generated. Thus with the 
same prescription as for the primary electrons, the simula- 
tion follows the evolution of secondaries as well. The simula- 
tion employs a computational box of 50 h^^ Mpc on a side, 
a grid of 512^ cells and 256^ dark matter particles. With 
the cosmological parameters described at the beginning of 
this Section, each numerical cell measures about 100 
kpc (comoving) and each dark matter particle corresponds 
to 2 X 10^ M0. Finally, momentum space for the CR 
components is subdivided into a few equidistant logarithmic 
momentum bins as detailed in Miniati (2002). 



3 THE 21CM RADIATION FROM THE 
INTERGALACTIC MEDIUM 

The 21 cm hyperfine transition of neutral hydrogen in the 
IGM provides a powerful probe for studies of the history of 

reionization. The physics of this transition has been well 
studied in the cosmological context, initially using linear 
theory (Scott & Rees 1990; Kumar et al. 1995; Madau, 
Meiksin & Rees 1997; Tozzi et al. 2000) and more recently 
using numerical simulations of reionization (CM; Furlanetto, 
Sokasian & Hernquist 2004). In this work we use the results 
from the simulations by CM (which we refer the reader to for 
further details) to compare the foreground to the intrinsic 
21 cm signal. 

The emission in the 21cm line is governed by the spin 
temperature, Tg. In the presence of the CMB radiation, Ts 
quickly reaches thermal equilibrium with the CMB temper- 
ature, TcMB, and a mechanism is required that decouples 
the two temperatures. While the spin-exchange collisions be- 
tween hydrogen atoms proceed at a rate that is too small 
for realistic IGM densities, Lya pumping contributes signif- 
icantly by mixing the hyperfine levels of neutral hydrogen 
in its ground state via intermediate transitions to the 2p 
state. When the first objects form, Lya pumping will effi- 
ciently decouple Ts from Tcmb if Ja > 10^^^ ergs cm~'^ s~^ 
Hz^^ sr^^. CM find that the expected diffuse fiux of Lya 
photons produced by the same sources responsible for the 
IGM reionization, satisfies the above requirement during the 
'grey age', from redshift ~ 20 to the time of complete reion- 
ization. As the IGM can be easily preheated by primordial 
sources of radiation (e.g. Madau, Meiksin & Rees 1997; Chen 
& Miraldar-Escude 2003), the universe will, most likely, be 
observable in 21cm emission at a level that is independent 
of the exact value of Ts. Variations in the density of neutral 
hydrogen (due to either inhomogeneities in the gas density 
or different ionized fraction) should appear as fiuctuations 
of the sky brightness of this transition (e.g, as a brightness 
temperature increment with respect to the CMB at an ob- 
served frequency !^ corresponding to a redshift 1 + z — Uo/f, 
given by 5T{v) = {Ts — Tcmb)t/{1 + z) where r is the opti- 
cal depth of a patch of IGM in the hyperfine transition; see 
e.g. Field 1959). In principle, high resolution observations of 
the 21cm line transition in both frequency (hence redshift) 
and angle should provide a detailed map of reionization. 

CM studied the evolution of brightness temperature 
fluctuations in both the 'late' and 'early' reionization mod- 
els (the S5 and the L20 models respectively, described in 
Section 2). The S5 model predicts a peak in the ampli- 
tude of the surface brightness fluctuations at ~ 115 MHz 
whereas in the L20 the peak occurs at f ~ 90 MHz. The 
overall amplitude of the signal at its peak, in both models 
is {ST^rnsY^'^ ~ 10 — 20 mK at angular scales S ~ 5 ar- 
cmin (see CM for details). When comparing the foreground 
fluctuations to that of the primary redshifted 21cm we will 
mainly use the S5 model where the peak in the 21cm signal 
occurs at frequencies which will most easily be probed by 
the upcoming telescopes (this choice will not affect any of 
the analysis performed here). We will examine the frequency 
and angular dependence of 21cm signal and foregrounds in 
the Conclusions. 
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4 EXTRAGALACTIC FOREGROUND 
SOURCES 

In this Section we describe the properties of the sources 
responsible for the extragalactic foreground signaL In par- 
ticular, we describe how we model their redshift distribution 
to calculate their spatial correlation functions and luminos- 
ity functions from the simulations. We will show that our 
models are able to reproduce reasonably well the observed 
populations and hence should provide us with a physically 
grounded framework for examining the contribution to the 
temperature fluctuations due to unresolved sources. 

We consider two main classes of sources, (a) Point 
sources: (i) radio galaxies/ AGN which produce synchrotron 
emission in their nuclei; (ii) HII regions inside ionizing halos 
which produce free-free emission, (b) Extended sources: (i) 
HII regions in the IGM which produce free-free emission; (ii) 
cluster radio halos and (iii) cluster radio relics, which both 
produce synchrotron emission. As we shall see in Section 5, 
radio galaxies, cluster radio halos and relics appear to pro- 
vide the most significant contribution to the sky brightness 
temperature fluctuations at angular scales of a few arcmin. 



4.1 Radio Galaxies 

In this Section we describe the simple scheme we have imple- 
mented in the simulations to reproduce the properties of the 
observed radio galaxy population. The physical processes 
that determine the amount of radio power from active galac- 
tic nuclei (AGN) are poorly understood. However, the radio 
power in AGN must be linked to the presence of the super- 
massive black hole in center of galaxies and to the accretion 
of gas onto it. Several groups have discussed the link between 
cosmological evolution of AGN and the formation history 
of galaxies using semianalytic models of galaxy formation 
(see e.g.; Monaco et al. 2000; Kauffmann & Haehnelt 2000; 
Cavaliere & Vittorini 2000; Wyithe & Loeb 2002; Volonteri, 
Haardt & Madau 2003; Haiman, Ciotti & Ostriker 2004 and 
references therein) and recently using cosmological hydrody- 
namic simulations (Di Matteo et al. 2003ab). Such models 
have been shown to successfully reproduce the B-band lu- 
minosity functions of QSOs and in some cases the X-ray 
luminosity functions as well. Here we make the hypothesis 
that the radio activity is also linked to the growth and ac- 
cretion of gas onto black holes and hence to the activity 
seen in the optical and X-ray bands. In particular, we use 
a prescription similar to that adopted in previous works to 
couple the cosmological simulations described in Section 2 
with a simple scheme for the radio population. Note that 
we shall not attempt to build a detailed model for the radio 
galaxy population but merely reproduce its main observed 
properties such as the luminosity functions and their angular 
clustering. 

The fraction of gas accreted by the central black hole in 
galaxies must be related to the observed correlation between 
stellar velocity dispersion and black hole mass (Ferrarese & 
Merritt 2000; Gebhardt et al. 2000). Using the most recently 
revised set of black hole mass and velocity dispersion mea- 
surements (Tremaine et al. 2002), the A/bh— o" relation gives: 

^- = (^-^±°-2)^^°'^-(2om^)'' 
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Figure 1. The radio luminosity function at 150 MHz in the sim- 
ulations (diamond symbols) plotted with Poisson error bars com- 
pared to the analytical fits (solid lines) of the observed luminos- 
ity functions from Willott et al. (2001) from the 3CRR, 6CE and 
7CRS samples. 

where a is the velocity dispersion of stars in the bulge. Fer- 
rarese (2002) and Baes et al. (2003) find that the stellar ve- 
locity dispersion measured in galaxies is strictly correlated 
with the asymptotic value of the circular velocity, Vc- 

log 14 = 0.9 log cr + 0.3. (2) 

KVc — (J relation consistent with the observations has also 
been measured from the simulations (Eq. 8; Di Matteo et 
al. 2003a). In order to avoid introducing any details on how 
gas cooling, star formation and feedback is linked to black 
hole fueling, we simply combine the previous two relations to 
find the mass accreted by a supermassive black hole. We as- 
sume this mass to be accreted at the Eddington rate after a 
timescale tq ~ tdyn ~ 0.1r„ir/Vc, where tdyn is the dynam- 
ical timescale estimated at a tenth of the virial radius, r^ir- 
All galaxies are assumed to undergo active phases with a 
duty cycle given simply by fq = tq/tH^z) where t„(^z) is the 
Hubble time at redshift z. This fraction gives us the expected 
number density of active black holes at a given redshift. We 
then use the empirical relation found by Merloni, Heinz & 
Di Matteo (2003) and Falcke, Koerding & MarkoflF (2003) to 
infer the radio luminosity. For each galaxy of a given nuclear 
X-ray luminosity and black hole mass, the relation gives, 

log L^"^ = 0.6 log Lx + 0.8 log Mbh + 7, (3) 

where Lx is the X-ray luminosity in the 2 — 10 keV band, 
which is related to the bolometric luminosity with a cor- 
rection Lx ~ O.OILboI (Elvis et al. 1994; Elvis, Risaliti & 
Zamorani 2002). In Eq. 3, is the radio luminosity at 

5 GHz. In Fig. 1, we show the radio luminosity functions at 
150 MHz and for z = 3, 2, 1,0.1 obtained from the simula- 
tions (diamond symbols) using in Eq. 3 and assuming 
a power law spectrum which gives us the luminosity of the 
radio sources at a given frequency simply as L^"^ oc v~" 
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Figure 2. Angular correlation function of the simulated radio 
galaxies (dotted line) compared to a power law with slope (in 
log-log) 1 — 7 = 0.8 and amplitude A = 2 X 10~^. The amplitude 
and shape of the angular correlation function is consistent with 
the one inferred from observations of faint radio galaxies in the 
NVSS and FIRST surveys. 



with spectral index a = 0.8. The results from the simu- 
lations are compared with the luminosity functions at 150 
MHz most recently measured by Willott et al. (2001) from 
the 3CRR, 6CE and 7CRS radio surveys (shown with the 
solid lines) § consistent with the models of Dunlop & Pea- 
cock (1990). Note that because of the limited size of our 
simulation box our luminosity function does include a small 
number of bright (very massive) objects. 

By projecting the maps from the simulations we also 
calculate the angular correlation function, w{6) of the radio 
galaxy population (shown in Fig. 2) defined as, 

wie) = {S{f)S{^+9)), (4) 

where 0, 6_ are angular coordinates, 9 — \6_\, and 5(0) = 
[n{(j)) / (n)] — 1, where n((;/)) is number density of galaxies per 
unit solid angle at angular position 0. Our estimate of w{9), 
is calculated from the Limber projection of the spatial cor- 
relation function of radio galaxies. The details of this calcu- 
lation are described in Section 5.4. The dash-dotted line in 
Fig. 2 shows the conventional power-law form w{d) = A9^''^ 
(e.g. Peebles 1993) where 1 - 7 = 0.8 and A = 2 x 10"^. 
The amplitude A and the shape of the angular correlation 
function are fully consistent with the values measured for 
cosmological clustering of the 1.4 GHz NVSS and FIRST 
radio surveys in the millijankys and sub-millijankys popu- 
lations which show that Aots is fairly constant at these fiux 
limits and Aots ~ 1 - 4 x 10"^ for 1 - 7 = 0.8 (Maglioc- 
chetti et al. 1998; Richards 2000; Overzier et al. 2003; Will- 
man et al. 2003). In summary, our simple model of the radio 
galaxy population broadly reproduces the observed luminos- 
ity functions and angular clustering of the faint radio source 
population. 



§ Note that the bumps seen in the fit to the luminosity function 
of Willott et al. (2001) are probably not real and may be due to 
restricted number of parameters in their model. The important 
feature is a break in the high and low luminosity populations with 
the low luminosity population disappearing at low redshifts. 



4.2 Free-Free Emission from HII Regions - ISM 
and IGM 

In this Section we will discuss the free-free emission from 
ionized regions, both in the ISM and IGM and how they are 
derived from the simulations. The free-free luminosity at 
frequency v from an ionized region of proper volume Vion is 
given by LI^ = t^Vion, where is the free- free emissivity for 
an hydrogen plasma at temperature T and electron number 
density rie, = 3.2 x lO-^^^n^ {T /IQ'^K)-^-^^ erg s'^ cm'^ 
Hz~^ (e.g. Rybicki & Lightman 1979; Oh 1999). 

The free- free luminosity associated with the ionized ISM 
of a galaxy can be estimated by taking the free-free emis- 
sivity to be directly proportional to the recombination rate. 
In particular, following the same prescription used by Oh 
(1999), we can write the recombination rate as: 

Nrec = aB-n^Vion ~ (1 — fe3c)Nion, (5) 

where as is the case B recombination coefficient (appro- 
priate for the densities of interest here), /eac is the escape 
fraction of ionizing photons and Nion is the ionization rate of 
the galaxy. Using the above equation to estimate the prod- 
uct n^Vion and assuming that the temperature of the ionized 
regions is T ~ 10^ K, it is straightforward to derive the free- 
free luminosity of each galaxy as: 

Li>''^'' = 3.2 X 10-39 (l^/-) jv.„„ erg s'^R^-'. (6) 

Q-B 

Once the galaxy emission properties and the value of /esc are 
set as described in Section 2, the value of Ll^''^'^ is easily 
determined for every object. As the mode of star formation 
is directly related to the gas metallicity, we expect a change 
from massive Pop III to more standard Pop II stars once 
the metallicity is Z > 10~'^Zq (e.g.; Bromm et al. 2001; 
Schneider et al. 2001). This would result in a lower efficiency 
for production of ionizing photons. To roughly take this into 
account we assume that for z < Zion, where Zion ~ 8 (13) for 
the S5 (L20) case, the spectrum switches to a more standard 
one, typical of Pop II stars. 

The emissivity from the more difi'use IGM HII regions 
can be directly derived from the reionization maps described 
in Section 2. From these, the value of Ue and the corre- 
sponding emissivity tl^'^'^'^ are known for each cell of the 
simulation box. After reionization, el^'"^^' can be derived 
assuming that the ionization fraction in each cell is unity. 

4.3 Galaxy Cluster Radio Halos and Relics 

Diffuse, nonthermal radio emission extending over Mpc 
scales is a well established observational feature for a fair 
fraction of galaxy clusters. Difi'use radio sources are classi- 
fied as either radio halos when their morphology is regular 
and typically centered on and resembling the X-ray emissiv- 
ity; or as radio relics when they are irregular and located 
at the periphery of the cluster (in this case the radiation is 
typically polarized). In general the volume integrated emis- 
sivity is well characterized by a steep spectrum, jiu) oc v~" 
with a > 1 (see e.g. Feretti 2002 for a review). 

Statistically it is found that roughly 30 — 40% of rich 
galaxy clusters (X-ray luminosity Lx > 10*^ erg s~^) do 
host a radio halo. In addition, it appears that the radio 
power (e.g. at 1.4 GHz) of radio halos scales steeply with the 
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Figure 3. Maps of the synchrotron emission (at 150 MHz and in 
units of erg cm~^ Hz~^ arcmin"'^) from shock accelerated 
CR electrons, i.e. so called radio relics at redshift 2=3 (top left), 
1 (top right), 0.5 (bottom left) and (bottom left). The map was 
obtained by adding the projections of the volume emissivity along 
the three coordinate axes. Thus, it is 50 Mpc on a side and 
150 Mpc deep. 



virial temperature of the associated clusters, i.e. Li,4ghz oc 
T'^ with /3 > 4 (Liang et al. 2000). The fraction drops below 
10% when considering smaller clusters (i.e. Lx < lO**^ erg 
s~^). It is generally believed, although still an open question, 
that this is mainly due to the limited sensitivity of current 
radio telescope (see, e.g., Liang et al. 2000; Bacchi et al. 
2003) 

The simulations described in Section 2 (Miniati et al. 
2001) successfully reproduce most of the observed properties 
of radio halos and radio relics, particularly those relevant for 
the purpose of this paper. Thus, diffuse radio emission pro- 
duced by synchrotron radiation from secondary e^ is char- 
acterized by a morphology and a low frequency radiation 
spectrum consistent with those of radio halos. In addition, 
the radio power (at 1.4GHz) scales with the cluster virial 
temperature with an index (5 ~ 4.2 which is consistent with 
the observations. Similarly, synchrotron emission from shock 
accelerated (primary) electrons is characterized by a spec- 
trum, morphology and location that resemble closely those 
of radio relics. 

For CR electrons with a distribution function N['y), 
such that N{'y)d'y is the number of particles with Lorentz 
factor 7 between 7 and 7 -I- d'y, the synchrotron emissivity 
(in units of erg cm"'' s~^ Hz~^) is computed according to 
the expression (e.g. Rybicki & Lightman 1979) 



F 



2^7" 



Svb sin a 



Ni-y) dy, 



(7) 



i^B = eB / {2-Kmsc) , is the electron nonrelativistic gyrofre- 
quency, F{x) = x dyK^/3{y) with A's/s the modified 
Bessel function of order 5/3, and all other symbols have 
their usual meaning. 

For the diffuse cluster radio halo emission, the popu- 
lation of CR electrons, N{'y), is provided by secondary e*. 
Consistently with current observational constraints, we con- 
sider that only 1/3 of galaxy clusters have radio halos, and 
further assume that radio halo formation is suppressed for 
objects of virial temperature below 1 keV. These two as- 
sumptions give us a fairly conservative estimate for this com- 
ponent. 

For cluster radio relics, on the other hand, the emit- 
ting electrons, N('y), are those directly accelerated at cosmic 
shocks. Even though observations are sparser for radio relics 
than radio halos, we find that the radio power of the simu- 
lated relics is consistent with the values inferred from obser- 
vations for a broad range of cluster temperatures (Miniati 
et al. 2001). In Fig. 3 we show four slices through the simu- 
lations of the emission from radio relics. We choose to show 
the maps for this component in particular to illustrate the 
complexity of its morphology. This will be important when 
discussing (in Sections 5.3 and 5.4) foreground removal. 



5 ANGULAR CLUSTERING OF 
FOREGROUND FLUCTUATIONS 

As we shall show in Section 5.4, the angular brightness tem- 
perature fluctuations due to foreground sources depend on 
two main quantities: the volume-average comoving emissiv- 
ity and their projected two point correlation function. We 
will start by discussing these quantities for each of the fore- 
ground sources described above. 



5.1 Comoving Emissivities 

For the extended sources and point sources respectively the 
volume average comoving emissivity is simply given by, 



e„(r, z)dV 



-y 

point 



- 17? LA^) 



(8) 



where B is the magnetic field strength, a is the parti- 
cles pitch angle (over which the ef,^" is further averaged). 



where V is the appropriate comoving simulation volume, 
ei,{r, z) in Eq. 8 is the emissivity at a given point in space, r, 
and Lv is the source luminosity which, for extended sources, 
is defined as = J e,ydVeff with the integral performed 
over the effective volume of the extended source (e.g. the 
cluster). The integration (summation) is performed over all 
points (sources) up to an emissivity (luminosity) correspond- 
ing to Li, < L,y{Sc,z) = Sc X 4,-Kd\{z) where Sc is the fiux 
cut off above which sources can be resolved and removed 
and dL the luminosity distance (cf. Section 5.3). 

In Fig. 4 we show the mean emissivity per comoving vol- 
ume at the observed frequency v — 115 MHz as a function of 
redshift for radio galaxies (dotted line), ISM free- free emis- 
sion (S5 run; long-short dashed line), IGM free- free emission 
(S5 run; short dashed line) , cluster radio relics (long dashed 
line) and cluster radio halos (dashed-dotted line) as obtained 
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Figure 4. Redshift evolution of the mean comoving cmissivity, 
at i/ = 115 MHz, due to: free-free from ISM HIl regions (long- 
short dashed lines); free-free from IGM HIl regions (short dashed 
lines) both in the S5 simulation; radio galaxies (thin dotted line 
and thick dotted line for Sc = od and Sc =0.1 mjy respectively); 
synchrotron from virialization shocks in radio relics (long dashed 
line); synchrotron from secondary electrons in radio halos gener- 
ated in cosmic rays (thin and thick dash-dotted line for Sc = oo 
and 5c = 0.1 mJy respectively). 



from the simulations. The thicker dotted and dash dotted 
lines show the emissivity when sources above Sc ~ 0.1 mJy 
have been removed. Here we will describe the total emissivi- 
ties from the various components and reserve the discussion 
of source removal to Section 5.3. 

As shown in Fig. 4, the radio galaxy, radio halos and 
relics were modeled in the simulation out to redshift z — 5, 
where most of the contribution is expected to come from. 
Below z = 5, the radio galaxy comoving emissivity domi- 
nates all the other components. Its peak, at z ~ 2, is at- 
tributed to the bright component of the luminosity function 
which is strongest around this redshift. However, the overall 
evolution of the radio galaxy {e^(z)) is fairly flat. In con- 
trast, for cluster radio halos and relics {e,^{z)) has a much 
steeper growth, by a factor ^ 10*, from z = 5 to z — 0. 
This is due to the rapid increase in the energy dissipated in 
shocks during merging/accretion events as a result of non- 
linear breaking of larger wave perturbations toward lower 
redshifts (Miniati et al. 2000). Finally, the free-free emis- 
sion from ionized regions is significant only at high redshift, 
when the other contributions become negligible. The emis- 
sivity from IGM HII regions increases at z < 20, as ion- 
ization proceeds. Once it is completed, at z ~ 8 for the S5 
case shown here, the emissivity decreases as the gas den- 
sity decreases. The emissivity from ISM HII regions is in- 
stead proportional to the star formation rate. The drop ob- 
served in the curve is related to the drop in the ionizing 
photon production rate following reionization, as described 
in Section 4.2. This emissivity is consistent with the estimate 




r (Mpc/^) 



Figure 5. The spatial correlation functions ^{r,z) for clus- 
ter radio relics, halos, radio galaxies (and ISM free-free emis- 
sion) and IGM free-free emission from top to bottom panel re- 
spectively. The spatial correlation functions are shown at z = 
0.1,0.5,1,2,3,5 (solid line, dotted, short dashed, long dashed, 
dot-short dashed, dot-long dashed) for radio halos and relics, 
at 2 = 0,1,3,5,10,13,19 (solid line, dotted, short dashed, long 
dashed, dot-short dashed, dot-long dashed) for the radio sources 
(and ISM free-free emission), and at z = 13, 14, 15, 16, 17, 18 for 
the free-free emission from the IGM (from bottom to top lines). 



of Bianchi, Cristiani & Kim (2001), which is derived inde- 
pendently based on the observed star formation rate and a 
population synthesis model. In any case, the contribution 
from ISM HII regions is always greater than that from IGM 
HII regions, with the ratio reaching a maximum of ^ 700 
at 2 ~ 20, when the IGM is mainly neutral, and a mini- 
mum of a few at lower redshift. In summary, we have shown 
that the radio galaxies, cluster radio relics and halos provide 
the largest contributions to the total comoving emissivity at 
moderate/low redshifts. The resulting brightness tempera- 
ture fluctuations however, are also a function of the spatial 
clustering of the various components (described in Section 
5.2). 
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Figure 6. The spatial correlation function for radio galaxies and 
radio halos at z = 2 and z = 1 respectively. The solid lines show 
the correlation functions of the total sample, the dashes those for 
sources below the flux limit Sc = 0.1 mjy. 



5.2 Spatial clustering 

We quantify the spatial clustering of the various sources 
via the redshift dependent two-point correlation function, 
^(r, 2), defined as (Peebles 1993): 



{(5(x,2;)(5(x + r, 2)), 



(9) 

£^(x,2)/(e^(2)) - 1 



where for extended sources, 5(x, 2) 
and ei/(x, 2) is the emissivity at a point with coordinate x at 
redshift 2. In practice, to evaluate ^(r, 2) we have computed 
the power spectrum Ps of (5(x, 2), and Fourier transformed: 



extended — 



(27r) 



„ ,, , sinfcr , , 2 ,, 

Ps{k) Ank dk. 

kr 



(10) 



For the point source populations (radio galaxies and free- 
free emission from ISM) 5(x, 2) = n(x, 2;)/(n(2)) — 1, where 
n(x, 2) is number density of galaxies at x. We evaluate the 
correlation function using the standard pair estimator. 



N{r,z) 



(11) 



^randomij' ^ 2) 

where N{r, z) is the number of pairs separated by a distance 
r, and Nrandom{r, z) is the number of pairs in the same r- 
bin expected if galaxies were randomly distributed. Because 
our simulated box is not periodic, we generate random cata- 
logues to compute Nrandomif, z) . Froui top to bottom, Fig. 5 
shows the computed ^(r, 2) for the cluster radio relics , £,rr, 
halos, £,RH, radio galaxies , (,rg, (the same as free- free emis- 
sion from ISM, ^ffjsM as they are the same halos) and 
free-free emission from IGM, ^ff.ioM (see caption of Fig. 5 
for details) . According to our model, the radio galaxies con- 
stitute a random sample of all halos (the duty cycle is not 
mass dependent, c.f. Section 4.1). As we expect, the corre- 
lation function of the galaxy population is close to a power 
law of the form ^(r) ~ {r/ro) ' (usually slightly suppressed 
on scales r ^ 20 — 50 kpc) and with a correlation length. 



ro which varies in the range of ro ~ 2 — 5/i~^Mpc in the 
redshift range studied (Fig. 5). Note that the amplitude of 
the correlation function grows with decreasing redshift for 
2 < 10. At higher redshift the trend is reversed probably as 
an effect of the mass resolution (when halos of mass corre- 
sponding to the resolution limit start representing the high 
end tail of the mass function). This is not important as all 
of the signal from point sources is dominated by relative low 
redshifts (c.f.; Section 5.4). 

The correlation functions for radio halos and radio relics 
have typically large amplitudes, a rather flat inner region 
and a steep decline at r > l/i~^Mpc (particularly in the 
case of radio halos). The correlation length for cluster radio 
halos is set predominantly by the size itself of the brightest 
sources, which is roughly of order of a Mpc or so. Radio 
relics are typically more extended, as reflected also by their 
correlation functions, which decline less steeply than in the 
case of radio halos. This is a consequence of the fact that 
accretion shocks can cover larger regions than that of the 
cluster, particularly when two or more groups/clusters are 
approaching before a merger (see Fig. 3). 

Finally, for the correlation functions of the IGM free- 
free emission we only show those corresponding to 2 > Zion- 
We note that as ionization proceeds the HII regions, which 
initially are isolated and surrounding the galaxies, will be- 
come increasingly more extended and, as expected, will start 
overlapping. This has the consequence that the correlation 
amplitude decreases with decreasing redshifts. Because the 
IGM free-free emissivity shown in Section 5.1 and its corre- 
lation function are negligible with respect to the other com- 
ponents (and in particular with respect to free-free emission 
from the ISM) we will not consider it further. 



5.3 Bright source removal 

As part of the actual experiment, bright enough sources will 
be directly identified as foreground contaminants, and hence 
removed from the observational maps without contributing 
to the brightness temperature fluctuations. The flux limit 
above which sources can be removed, Sc, will of course be 
dependent on the instrument sensitivity. Here we take Sc = 
0.1 mJy, corresponding to the expected sensitivity for source 
detection by LOFAR and discuss the major effects that this 
has on both the emissivities and correlations functions. 

Note that removal of point sources will be far easier than 
the removal of extended components. Here we shall illustrate 
the effects of applying a flux cut on radio galaxies and on 
cluster radio halos (because, as we shall show in Section 5.4, 
these are the two dominant contaminants of each class). For 
cluster radio halos, we go back to the volume emissivity dis- 
tribution calculated from the simulation and calculate the 
total radio luminosity of each identifled cluster by integrat- 
ing the synchrotron emissivity over its volume (roughly out 
to the virial radius, r„ir), Lv = Jq^"^^ ei^Aivr^dr. The cluster 
radio emission is then removed from the volume emissivity 
distribution if > Sc x A-Kd\. As shown by Fig. 3 and dis- 
cussed in Section 4.3 the morphology of radio relics is fairly 
complex. We will not attempt here to remove this extended 
component. However, we do expect it to provide similar ad- 
vantages as those gained from the removal of cluster radio 
halos which we describe below. Nevertheless, such an oper- 
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Figure 7. Prediction for the correlation signal owing to intensity 
fluctuations of radio galaxies (dotted lines), ISM free- free emission 
(short dashed), radio relics (dashed dotted), radio halos (long 
dashed) at = 115 MHz. The thicker lines show the signal when 
sources above a flux 5c = 0.1 mjy are removed. The solid line 
shows the primary correlation signal due to the redshifted 21cm 
emission (CM) 



For the radio galaxies, we see from Fig. 6 that the corre- 
lation function of objects below a flux cut is biased in a scale 
dependent way with respect to the whole population. This 
is very different to the alternate case of galaxies brighter 
than flux Sc- For the latter, we find that the correlation 
function of the brighest galaxies is similar in shape to that 
of the whole population (linear biasing) , even on scales sig- 
nificantly smaller than the correlation length (as found by 
e.g.. Mo & White 1996). Note that we have checked that the 
effects we have just described, scale dependent bias versus 
linear bias, are the same when a mass cut is adopted in- 
stead of a ffux cut. We have checked this by calculating the 
correlation function ^ for objects below and above a given 
dark mass limit Mmin- This is expected as in our model 
we have directly linked e.g., the radio galaxy population to 
the dark matter potential. The relative change in clustering 
of halos after applying an upper mass threshold is not often 
considered, but the boost in clustering on small scales which 
results, along with suppression on linear scales, is expected 
in hierarchical models of structure formation (e.g. Sheth & 
Tormen 1999). 

Finally, we note that, as we shall show in the next Sec- 
tion, the dependence of the spatial clustering on the ffux 
cutoff is important when the spatial functions are mapped 
into the angular clustering signal. Such an effect had not 
been take into account previously in calculations of sky in- 
tensity fluctuations due to clustering of foreground sources 
within the context of the 21cm experiment. 



ation for cluster radio relics may turn out more difficult due 
to their more irregular morphology. 

As shown in Fig. 4, the removal of sources above the 
detection threshold has the obvious effect of decreasing the 
volume emissivity (c.f; Eq. 8). Because of the relatively low 
ffux cut level the emissivity drops dramatically towards the 
low redshifts. In particular, for radio galaxies emissivity de- 
creases between a factor of ~ 10 to ~ 200 from « = 5 to 
2 = 0.1 with respect to the case when no sources have 
been removed. For radio halos removal of sources above 
Sc = 0.1 mJy significantly affects the emissivity below z ~ 2. 
At z < 2 the emissivity decreases by up to a factor of a few 
thousand. 

In addition to affecting the overall emissivity, the re- 
moval of bright sources has an important effect on the spatial 
correlation functions. As an illustrative example, we show in 
Fig. 6 the correlation function for radio galaxies and radio 
halos at z = 2 and z = 1 respectively (solid lines) com- 
pared with the case in which a ffux cut, Sc = 0.1 mJy is ap- 
plied (dashed lines), and all the objects with ffux exceeding 
Sc have been removed from the respective simulated cata- 
logues. Typically, removing the bright sources has the effect 
of steepening the spatial correlation function and hence de- 
creasing the correlation lengths. In particular, for the radio 
galaxies, the amplitude at scales r ^ l/i~^Mpc is boosted, 
but at scales larger than the typical correlation length it 
is significantly suppressed. For the case of radio halos, the 
effect is similar and the suppression of the correlation func- 
tions occurs at even smaller scales then in radio galaxies. 
Also, because the correlation functions of these sources de- 
cline very steeply the overall effect is more dramatic at large 
scales. 



5.4 Angular Clustering 

We now calculate the contribution to the intensity ffuctua- 
tions, 51, owing to the inhomogeneities in the space distri- 
bution of the various sources described above. The angular 
correlation function of such intensity ffuctuations is given 
by, 



C(e) = {5I{^)SI{<f> + 9)), 



(12) 



where (j), 6 are angular coordinates, and 51{(j)) — [I{<l))/ (I)] — 
1, where /(</>) is the intensity at angular position (f). 

Following the standard approach, we use the Limber 
projection (Limber 1953; Peebles 1993) to calculate, C{6) 
(in units of erg^ s~^ cm"* Hz"^ sr~^): 

1 



c(e) 



(47r)2 



dVidV2 , , — r , ,2 , — r5(r, 2i,22),(13) 



where dVi,dV2, {eu{zi)), {e^{z2)) are the comoving volumes 
and emissivities at redshifts zi and 22 respectively. For small 
angles 6 and small radial distances, 21 ~ 22 ~ z and we 
write: 



1 



(47r) 



dz 



{€v{z))^ dx 

(l + z)2 



dui{r,z), (14) 



where ^(r, z) is the spatial two point correlation function, 
— + x^O^ , and x is the comoving coordinate, related 
to the redshift (for the ffat models) by: 



x=i- I [Vlo{l + z'f + VlA\''^''dz\ 
no 



(15) 



and we have written dz\ dz2 = (dz/dx) du dz for x\ ~ 
a:: 2 ~ a:: in Eq. 13. Analogously, using the standard Limber 
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Figure 8. Contours of the angular power spectrum (in mK) of the the 21cm fluctuations (shaded area) and the angular power spectra 
due to intensity fluctuations of radio galaxies (solid contours) and radio halos (dot-dashed contours) as a function of angle (/) and 
observed frequency. The lines in the right panel show the contours for the total foreground signals. In the left panel, foreground sources 
below a flux limit Sc ^0.1 mjy have been removed. 



equation the angular correlation function, 1^(9) (c.f., Eq. 4), 
can be written as, 



dz 

du dz — ^(r, z) 
dx 



(16) 



We can then simply rewrite the angular clustering of in- 
tensity fluctuations in terms of a brightness temperature 
fluctuation in the sky as: 



/y2 \ 

\-'-rms / 
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dT 



1/2 



'l{l + l)Ci' 
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\ dT J 



2k /kT 
c2 VIT 



(e^ 



(17) 



(18) 



with X = hiy/kT and T = 2.725 K. The power spectrum, Ci 
(in units of erg^ s~'^ cm~* Hz~^ sr~^), is related to C{0) by 
the Legendre transform as: 

Ci = 27r J d{cose)Pi{cose)c{e)^2-K deejoiw)c{9),{i9) 

where Jo is the Bessel function and the last step is valid for 
1^1. Fig. 7 shows the prediction for the correlation sig- 
nal caused by the intensity fluctuation due to radio galaxies 
(dotted lines), ISM free- free emission (short dashed), clus- 
ter radio relics (dashed dotted), cluster radio halos (long 
dashed) at i/ = 115 MHz, the frequency at which the 21cm 
emission peaks, according to S5 model of CM. The thinner 



lines show the total signal when no flux cutoff has been ap- 
plied, while the thicker lines show the signal when sources 
above a flux Sc — 0.1 mJy (corresponding to the planned de- 
tection limit of LOFAR) are removed. The solid line shows 
the primary correlation signal due to the redshifted 21cm 
emission. The total Trms signal due to any of the foreground 
components exceeds the primary 21cm signal. In particu- 
lar, the extended sources such as the radio halos and ra- 
dio relics provide fluctuations about three (and more) or- 
ders of magnitude and the radio galaxies about two orders 
of magnitude larger than the 21cm emission through the 
whole range of angular scales. The rms fluctuations due to 
free-free emission from the ISM, instead, are close to those 
of the 21cm. However, the effective brightness temperature 
fluctuations due to foreground sources after the removal of 
sources with S > Sc (for illustration we have applied the 
removal only to the highest point and extended components 
of the rms fluctuations, e.g. radio halos, and radio galax- 
ies), decrease signiflcantly. For the radio galaxies, detection 
and removal of sources above Sc ~ 0.1 mJy decreases the 
fluctuations due to this component down to the level of the 
21cm emission at scales 9^1 arcmin. Above these angu- 
lar scales the foreground signal drops out completely. This 
is a result of the combined effects of the decrease of the 
overall emissivity (e.g.. Fig. 4) and the suppression of the 
the amplitude of spatial correlation of faint radio sources 
at the scales that map into angular scales of 9 ^ 1 arcmin 
(where radial scales and angular scales are roughly related 
as, r ~3(6l/l')[(l + 2)/l]°-^/i"^ Mpc). For this flux cut, an- 
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gular scales of 9 ^ 1 arcmin are thus free of this contaminat- 
ing signal, rendering the 21cm tomography free of the radio 
galaxy contamination. A similar result is obtained when the 
flux cut is applied to the radio halos. The decrease in ampli- 
tude in this rms fluctuation is even more pronounced than 
in the case of radio galaxies. The signal drops to the level of 
the primary 21cm at ^ < 1 arcmin and is also well below it 
at larger angular scales. Applying the same flux cut to the 
ISM free-free emission sources, would lead to a similar efi^ect 
as for the radio galaxies, but in this case, the signal would 
drop well below that of the 21cm on all scales (and below 
the axis range plotted in Fig 7). 

After removing sources with flux above 5'c = 0.1 mjy, 
the doniiuant contribution to the rms fluctuations from radio 
halos comes from 0.6 < z < 2. Similarly for the case of 
the radio sources (which however show a different evolution 
of their mean emissivity and different clustering), the main 
contribution to the rms fluctuations is in the range 0.5 < 
z < 2.5. This is set by the the projection of the emissivity 
and correlation function as shown in Eq. 14. 



6 CONCLUSIONS 

We have examined the contribution of cxtragalactic fore- 
grounds to the angular brightness temperature fluctuations 
at the frequencies relevant for observations of the redshifted 
21cm signal. We have used a set of cosmological simulations 
to model the evolution of radio galaxies, free-free emission 
from the ISM and IGM and radio halos and relics. We cal- 
culated their comoving space density emissivity and spatial 
correlation functions and projected the simulation volumes 
to predict their expected sky temperature fluctuations. 

Because high resolution observations of the 21cm signal 
will be carried out in both frequency and angle it is best to 
summarize the main conclusion from our analysis by show- 
ing the intrinsic and foreground signals as a function of both 
scale and frequency. Fig. 8 shows the angular power spec- 
trum of the 21cm signal (from CM simulations), and that 
of the radio gala^xies and radio halos (the two dominating 
foreground components for which source removal was stud- 
ied here) as a function of frequency (where we have taken 
the spectral index of radio galaxies and radio halos to be 
a = 0.8 and 1 respectively) and angle (scale I). In summary, 

• As previously emphasized by Di Matteo et al. (2002, 
for the case of radio galaxies) and Oh & Mack (2002, for the 
case of free-free emission from the ISM) we have found that 
the absolute foreground signal fluctuates more strongly than 
that of the 21 cm signal at all angular scales and frequen- 
cies (left panel Fig. 8). In addition, here we have shown that 
the largest fraction of these fluctuations is due to extended 
cluster radio halos and relics, both of which had previously 
boon neglected in the context of the 21cm tomography ex- 
periment (see however Waxman & Locb 2000). Radio emis- 
sion from the cores of galaxies contributes more strongly 
than the free-free emission from their ISM. The presence of 
foreground signal many orders of magnitude above the pri- 
mary signal needs to be addressed because it may seriously 
hinder the detection of the 21 cm line emission as the lat- 
ter could be easily mimicked by small errors/irregularities 
in the substraction of the foreground spectra. 



• We have considered source removal above a flux limit 
Sc = 0.1 mJy, (close to LOFAR sensitivity) for point sources 
and radio halos (as an example of an extended source). We 
have found that, as expected, the amplitude of the fore- 
ground power spectra of both radio halo and radio galaxies 
decreases at all scales. In particular, with the additional spa- 
tial information on the foreground sources provided by the 
simulations, we have shown that the amplitude of the an- 
gular foreground fluctuations drops well below that of the 
21cm at scales 6^1 arcmin and becomes comparable to 
it at scales smaller that 9^1 arcmin (right panel Fig. 8). 
At angular scales larger than ^ ~ 1 arcmin and at frequen- 
cies v < 150 MHz (simply corresponding to the rcdshifts 
prior to full rcionizatioii in the S5 model) the 21cm signal 
should be detected directly. At scales smaller than 9 ~ 1 
arcmin, substraction of the foregrounds can be feasibly car- 
ried out in frequency space, as the overall foreground signal 
is at most of the same order as the primary 21cm. Com- 
paring maps closely spaced in frequency, (as discussed in 
detail by Zaldarriaga et al. 2004) after the removal of bright 
sources should make the detection of the signal feasible at 
small scales also (where the foreground power spectra are 
different than the one of the 21 cm; Fig. 8)^. 

From Fig. 8, wc conclude that bright sources contribute 
to the angular power spectrum at ^ > 1 arcmin and, once re- 
moved, the power on these scales is suppressed. This effect is 
due to a scale dependent bias of the spatial correlation func- 
tions of sources below a certain flux cut. This result empha- 
sizes the importance, when determining the angular fluctua- 
tion and power spectra of foregrounds, of using simulations 
to project the detailed spatial distribution of foregrounds 
as a function of flux cut. Analytical estimates of brightness 
temperature fluctuations (Di Matteo et al. 2001; Oh & Mack 
2002), which had to assume a certain angular clustering of 
the radio foregrounds, appear to be insufficient for determin- 
ing the effect of these contaminants on all scales. This had 
lead previous authors to generally more pessimistic conclu- 
sions for the feasibility of 21cm tomography. We flnd instead 
that, due to the effects of bright source removal, 21cm ob- 
servations carried out at angular scales 9 of the order of a 
few arcmin should not be strongly affected by cxtragalac- 
tic foregrounds. However, how source substraction can be 
carried out from the observed maps, needs detailed work. 
The maps we have simulated will need to be folded with the 
instrument response functions as a function of angle and fre- 
quency (and other specific details of the instrument design). 
This is currently being investigated for the case of LOFAR 
by Valdes et al. (2004). 

We note that determining the angular scale at which 
the foreground signal should drop out is, of course, sub- 
ject to some uncertainty associated to specifics of source 
removal. Besides needing to account for the instrument re- 
sponse, which is beyond the scope of this paper, this will 
of course be a function of the theoretical approach we have 

^ Wc note further that the amplitude of the 21cm signal shown 
in Fig. 8 corresponds to a given bandwidth Au = uAz/{l + z) = 
IMHz and is expected to increase with bandwidth resolution. 
CM show that the signal can increase by a factor up to 4 at small 
angles and high frequencies if e.g.; Au = O.lMHz can be achieved 
making foreground substraxition less crucial at scales 6^1' 
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adopted to model the foreground source populations down 
to flux limits currently unobserved. In fact, it is part of the 
scientific goals of future high sensitivity low frequency radio 
facilities such as LOFAR, PAST and SKA to sample for the 
first time faint radio galaxies and diffuse cluster radio omis- 
sion (radio halos and radio relics, see Enfilin & Rottgering 
2002, Briiggen et al. 2004) and hence determine the proper- 
ties of these population up to high redshifts. 
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